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Local  projection  methods  which  yield  C^m”^  piecewise  polynomials  of 
order  m+k  as  approximate  solutions  of  a  boundary  value  problem  for  an  m 
order  ordinary  differential  equation  are  determined  by  the  k  linear 
functionals  at  which  the  residual  error  in  each  partition  interval  is  required 
to  vanish.  We  develop  a  condition  on  these  k  functionals  which  implies 
breakpoint  superconvergence  (of  derivatives  of  order  less  than  m)  for  the 
approximating  piecewise  polynomials.  The  same  order  of  super-convergence  is 
associated  with  eigenvalue  problems.  *xA  discrete  connection  between  two 
particular  projectors  yielding  0(|A|  )  superconvergence ,  namely  (a) 

collocation  at  the  k  Gauss-Legendre  points  in  each  partition  interval  and 
(b)  "essential  least-squares"  (i.e.,  local  moment  methods),  is  made  by  askina 
that  this  same  order  of  superconvergence  result  when  using  collocation  at 
k-r  points  per  interval  and  simultaneous  local  orthogonality  of  the  residual 
to  polynomials  of  order  r;  the  k-r  points  then  necessarily  form  a  subset 
of  the  k  Gauss-Legendre  points. 
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SIGNIFICANCE  AND  EXPLANATION 


In  piecewise  polynomial  collocation,  one  approximates  the  solution  of  an 
ordinary  differential  equation  (ODE)  by  the  piecewise  polynomial  function  (on 
a  given  subdivision  of  the  interval  of  interest  and  of  a  given  degree)  which 
satisfies  the  ODE  at  a  certain  set  of  points,  the  collocation  points  (and 
satisfies  the  side  conditions).  In  1971,  the  authors  discovered  that  a 
careful  placing  of  the  collocation  points  (viz.  the  choice  of  the  appropriate 
Gauss-Legendre  points,  familiar  from  Gauss  quadrature,  in  each  interval  of  the 
subdivision)  would  achieve  an  unusually  high  order  of  convergence  at  the 
break-points  of  the  chosen  subdivision.  This  phenomenon,  which  was 
subsequently  dubbed  (by  Douglas  and  Dupont)  "superconvergence",  has 
contributed  much  to  the  present  popularity  of  collocation  methods. 

In  this  report,  the  authors  show  that  such  superconvergence  is  not  a 
special  feature  of  collocation  alone,  but  can  be  achieved  by  a  whole  class  of 
local  piecewise  polynomial  approximation  methods,  with  collocation  at  Gauss 
points  at  one  extreme  and  "essential  least  squares"  at  the  other.  In  this 
latter  method,  the  residual  error  is  not  made  to  vanish  at  certain  points  (as 
in  collocation),  but  is  made  to  be  orthogonal  to  a  space  of  piecewise 
polynomial  functions  (of  appropriate  degree)  on  the  same  subdivision. 

As  part  of  the  discussion,  the  authors  provide  an  argument  which,  so  they 
hope,  gives  more  insight  into  collocation  at  Gauss  points  than  the  original 
one. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MPG,  and  not  with  the  authors  of  this  report. 


LOCAL  PIECEWISE  POLYNOMIAL  PROJECTION  METHODS  FOR  AN  ODE 
WHICH  GIVE  HIGH-ORDER  CONVERGENCE  AT  KNOTS 


Carl  de  Boor  and  Blair  Swartz 

Introduction.  This  is  the  last  in  a  triple  (see  [2],  [3])  of  papers  concerned  with  high- 
order  approximation  to  eigenvalues  of  an  ODE  using  collocation  at  Gauss  points. 
Correspondingly,  its  two  sections  are  labelled  9  and  10  ,  but  it  can  be  read  without 
reference  to  [3],  i.e.,  to  Sections  5  -  8  .  Items  labelled  x.y  or  (x.y)  are  to  be  found 
in  Section  x,  e.g.,  in  [2]  in  case  x  is  less  than  5  . 

When  writing  [2] ,  we  were  forced  to  go  through  the  arguments  in  [1]  once  again  and 
ended  up  improving  upon  them  somewhat  (see  the  proof  of  Theorem  9.2  below).  In  the  process, 
we  considered  more  general  local  piecewise  polynomial  projection  methods  in  an  effort  to 
discover  just  what  produces  the  superconvergence  at  breakpoints  in  Gauss-point  collocation. 
This  led  us  to  a  simple  set  of  conditions  on  the  local  projector  used  which,  so  we  found, 
had  been  formulated  much  earlier  by  Pruess  (4)  in  another  context.  In  addition  to  updating 
our  earlier  results  in  [1]  and  (2]  to  cover  this  wider  class  of  projection  methods,  we  give 
a  detailed  analysis  of  these  special  local  projectors  and  establish  a  simple  link  between 
the  two  best  known  among  these,  viz.  Interpolation  at  Gauss  points  and  Least-squares 
approximation. 

9.  Sane  projectors  which  yield  super  convergence.  As  de  Boor  and  Swartz  [1]  describe  it,  / 
local  projection  methods  which  involve  sufficiently  rough  piecewise  polynomials  are 


basically  determined  by  a  bounded  linear  projector  2  which  carries  C [  —  1 , 1 ]  onto 
(polynomials  of  order  k  ,  i.e.,  of  degree  <  k),  hence  satisfies 
(1)  #f  -  Qfl  <  const  #Dkf II  ,  all  feC(l°[-1,1] 

V  00 

l 

for  some  constant  constg  .  Then,  given  a  partition  A  :=  (t  )Q  of  [0,1]  with 

0  “  tn  <  ...  <t.  =  1,  | A |  s  =  max  At.  ,  I 


I  s?''-  : 


Q  determines  a  map  projecting  ^x^ctt^  1»ti!  =!  onto  ^  (the  space  of 


Sponsored  by  the  United  States  Army  under  Contract  Nos.  DAAG29-75-C-0024  and  HAAG^-SP- 
C-0041  and  by  the  United  States  Department  of  Energy  under  Contract  No.  W-7405-Ena.  3r'. 


piecewise  polynomials  of  order  k  with  breakpoints  in  A  )  by  translating  the  procedure 
for  Ct-1,1]  to  each  partition  interval;  i.e.,  by  requiring  that,  on  each  an<*  f°r 


y  e  c4  , 

(2)  QAy  =  S~1QSiy  ,  with  (Sig)(s)  :=  g( t.+1 /2+sAt . /2 )  for  s  e  [-1,1]  • 

Then,  from  ( 1 ) , 

<V  ,f-«Af,<i>  <  const2|At.|kIIDkfll(i) 

with 

"gl(i)  :=  SU£>  !  ti<x<ti+1^  ' 

Finally,  the  projection  method  for  the  m-th  order  differential  equation  Mx  =  y,  Bx  =  0 

is  determined  by  requiring  that  xA  e  pm  :=  P  ,,  A  n  C^m  1 ^  [0,1]  satisfy 

u  m>  k  ,u  n\>  K  ,u 

(3)  V*A  =  e^y,  |x4  =  0  . 

We  consider  a  set  of  constraints  upon  Q  which  permits  proof  of  0(|A|k  nJ 
breakpoint  superconvergence  for  this  projection  method.  These  constraints,  constructed  by 
Pruess  in  another  context  [4,  pp. 553-4,  esp.  p.554,  line  5],  can  be  stated  as  follows: 


For  some  positive  integer  n  <  k  (and  in  terms  of  l^f-l/l] 

<4)  *i  1  (1-°)[Vn+1-i]  '  i=1 . "* 


This  condition  is  equivalent  to  the  following:  For  some  sequence  (f^)k  n  with 

j 


(5a) 

we  have 
(5b) 
and 
(5c) 


P.  =  span  ( f i ^  ,  all  j. 


Qfj  =  0  for  j>k 


/  V 


i  3 


for  i  <  k  <  j  4  k+n+1  -  i 


Indeed,  by  (5a),  (4)  is  equivalent  to  having 

[f  =  0  for  r  4  i  ,  j  <  k  +  n+1  -  i  and  i = 1 , . 

In  fact,  since  (1-2)f;  -  0  for  j  <  k  ,  f 4 )  is  equivalent  to  having 


/fr(1-Q)f 


0  for  r  <  i,  k  <  j  <  k+n+1  -  i  and  i 


-2- 


j 

i 
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i.e.,  for  j  >  k  and  r  <  k+n+1  -  j  ,  and  so,  by  (5b),  (4)  is  equivalent  to  (5c).  This 
shows  that  (5a-c)  implies  (4).  On  the  other  hand,  for  any  linear  projector  Q  onto  P^  , 
we  can  find  (f^)^+n  satisfying  (5a-b)  by  taking 

g.  ,  i  <  k 

f.  s*  1  with  (g.)  s.t.  P.  =  span(g  )3,  all  j  , 

(1-Q)gi  ,  i  >  k  1  3 

hence  the  argument  also  shows  that  (4)  implies  (5a-c).  Finally,  this  last  statement  shows 

(with  gA(t)  =  ti_1,  all  i)  that  (4)  is  also  equivalent  to 

(6)  tr(1-Q)tS  =  0  for  r  <  n,  r+s  <  k+n 

Since  Pruess  was  the  first  to  consider  projectors  satisfying  (4)  (or  (6),  see  [5)),  we 
call  any  linejr  projector  Q  onto  Pk  and  satisfying  (4)  a  ssuper  projector  of  order 
(k,n). 

Example  1.  Collocation  Taking,  in  particular,  f  (t)  »  H.tt-p.)  ,  with  p  ,  ...,  p 

1  J'l  ]  I  K 

the  collocation  pattern  and  p  ,  ....  p  arbitrary,  we  find,  from  (5),  the  condition 

k  + 1  k+n 

1  k 

/  p(t)  II  (t-p.)  dt  «  0  for  all  p  e  P 

-1  j=1  3  n 

(used  in  [1])  to  imply  that  <2  ,  given  by  polynomial  interpolation  at  p,  p,  isa 

ssuper  projector  of  order  (k,n). 

Example  2.  Essential  least  squares  (method  of  moments,  or  of  iterated  integrals) 

Taking,  in  particular,  f^  =  P^_i  *m  the  Legendre  polynomial  of  degree  i-1  ,  all  i  , 

we  find  that  Q  ,  given  as  Least  squares  approximation  from  P^  ,  is  a  ssuper  projector,  of 

order  (k,k)  .  We  have  called  the  corresponding  process  "essential  least  squares"  because 

the  associated  projection  method  (3)  requires  that  the  residual  error,  Mx^  -  y  ,  be 

orthogonal  to  P  =  d"1  [P™  ,  O  ker  S  )  (assuming  §  “  (P.)T  to  be  linearly  independent 

kpA  ro+kfA  —  —  l  i 

on  P^);  while  ordinary  least  squares  asks  that  this  residual  be  orthogonal  to 

M[pm  n  ker  S  1  •  This  process  has  also  been  called  a  "local  moment  method"  for  an  mth 
m+k ~ 

order  equation.  In  this  connection,  recall  that  Wittenbrink  [6;Ex.3c]  shows  this  to  be 
equivalent  to  asking  that  the  iterated  integrals  of  order  j  ,  1  <  j  <  k  ,  of  the  residual 
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error  vanish  at  all  the  breakpoints*  We  have  chosen,  however,  to  emphasize  its  connect! 
to  least  squares. 

The  validity  of  (4)  suffices  for  proof  of  the  following  result  from  which  we  shall 
k*tn 

conclude  C(|i|  )  breakpoint  superconvergenee.  The  lemma  (and  its  proof)  are  a  varia 

of  Pruess*  result  [4;  section  3]  and  [5;  Lemma  2] . 


Len»a  9.1  Let  Q  be  a  ssuper  projector  of  order  (k,n)  .  Then 

1 

<7>  I  /  f(1-Q)g  |  <  const  l  IID^f llDk+n"^gl  • 

-1  2  j<n 

Proof.  Let  (T.f)<t)  :=  E^.ohtO)  tX/i!  .  Then 

1  1 

/  f(1-Q)g  =  /  f  ( 1  -Q  )T  g  +  0  (  Ilf  II  IIDk+ngll)  , 

-1  -1  *  n 

while 


1  k+n-1  l 

/  f(1-Q)Tk+ng  =  I  Drg(0)/r!  /  f(t)(l-Q)tr  , 
-1  r=k  -1 


since  Q  reproduces  Pk  •  On  the  other  hand,  since  Q  is  a  ssuper  projector  of  order 
(k,n)  , 


-1 


(using  (6)),  and  so 

1  1 

/  f(1-Q)tr  =  /  (  f  -  Tk+n_rf  )  (1-Q)tr  =  D(  llDk+n_rf  II;  . 


-1 


Consequently, 


1  k+n 

/  f(1-Q)g  =  0{  l  IIDrgllllDk  +  n‘rfllJ 

~1  r=k 

and  the  substitution  j  :=  k+n-r  brings  this  into  the  form  (7).  ||] 


With  the  definition  (2)  of  Q  ,  it  follows  that 

1  £-1  1 
1/  f(i-s.)g  I  =  I  r  (dt./2>  /  (s.f)d-o)(s  c)  i 
o  A  i=o  1  -i  1  1 

£-1 

<  const  E  (&t./2)  Z  II oh  fHr>k 
Q  .  n  i  ..  i 

i=0 
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IID^fl  =  (Ati/2)jllD3fli{i) 


Consequently,  we  have 


Corollary  1  If  Q  is  a  ssuper  projector  of  order  (k.n)  ,  then  there  exists  const  _so 
l-1  2-1  ..  . 

that  for  f  e  x  ^ti<ti+1l  and  g  e  x  Lw  tti#ti+1]  ' 

0  0 

1/  f d-CA>g  I  <  const  I  (At/^Ufl  Igll 
0  1=0 


If  1  ...  :=  max  «DJf  II,  .  .  . 


If  now  f  and/or  g  in  Lemma  9.1  are  not  as  smooth  as  required,  say, 

(n  )  (k+n  ) 

f  s  L^  ,  g  e  9  ,  with  n^ ,  n^  <  n  ,  then  we  are  only  entitled  to  consider 

T,  g  and  T.f  for  j  <  n.  ,  hence,  instead  of  (7),  we  get 
K+n„  3  .  f  ..... 


1/  fO-0)g  )  <  const  l  llDjfllllD 

-1  8  j<n{ 


k+(n  -j>+ 

■  II  n  9 


Correspondingly,  we  get 


Corollary  2  Lf  Q  is  a  ssuper  projector  of  order  (k,n),  then  there  exists  const  £o 
2-1  <n  )  2-1  (k+n  .) 

that  for  f  £  xL.  '  It.  ,t.  and  qe  x  L_  9,1  It.,t.,,]  with  n.  .  ,  n  .  <  n  , 

— - - - —  0  i  1+1  -  0  ii+1  -  f,i  9,i  — 


all  i, 

J  2-1  1+k+min{n  . ,n  .} 

(9)  1/  f < 1 -Q  )g  I  <  const  I  (At  )  9,1 

0  i=0 


n f "  ,.,igiL  ...  . 

n_  . , ( l)  k+n  , , ( i) 
f,i  g, i 


We  now  sketch  proofs  concerning  the  convergence  of  x^  satisfying  (3)  to  x  .  From 

the  proof  of  Theorem  3.1  in  de  Boor  and  Swartz  [1],  we  find  that  x^  exists  uniquely  for 

|A|  sufficiently  small;  that  if  Mx  E  C  [0,1]  ,  then 

BDr(x-x,)ll  <  const  u>.  (Mx),  0  <  r  <  m  , 

with  w^tf)  •-  supj  sup(  |  f  ( t)-f  (  s)  |  :  t^<s,t<t^+1)  ;  and  that  if 

x  E  C(m  11  [0,11  C^k)  [0,1]  ,  then 

R  D r  ( x— x , )  II  <  const  |A|k#xll  ,  0<r<m  ;  llxll  .  ,  :=  max.  8x1.  . 

A  “  m+k , A  ] ,A  l  j ,  (  i ) 

The  proof  of  Lemma  4.1  in  that  paper,  which  uses  this  convenience  of  x.  =:  Rx  together 


*  fiaai  (Hr ilr f  Ma  n  i*  i 


|Dr(x  -  x  )(t.)|  <  const  |A|k+n!lxll  ,  0  <  r  <  m  . 

A  1  m+k+n,A 

Remarks.  Isolated  solutions  in  m+k+n^ [0 r 1 ]  to  nonlinear  problems  can  be  handled  as 

in  [1;  Thm.  3.1],  where  the  question  of  superconvergence  is  reduced  to  the  superconvergence 

associated  with  a  linearized  problem  (which  we  have  just  settled) . 

We  have  left  open  so  far  the  question  of  when  the  side  conditions  0  are  "suitable". 

Simply  put,  the  side  conditions  are  "suitable"  if  Green* s  function  resulting  from  it  allows 

the  earlier  argument  to  be  made.  If,  for  example,  0  consists  of  multipoint  conditions, 

then  one  fixes  a  partition  A^  =  of  [0,1]  whose  partition  points  contain  all  the 

points  involved  in  0  ,  and  insists  that  all  partitions  A  under  consideration  are 

refinements  of  A^  •  Green’s  function  for  (M,  0)  then  satisfies 

(10)  (3G/3sr)G(s,-!  e  c'n)(0,S]xc'n)[s,1!  , 

0  0 

and  this  is  enough  to  complete  the  argument  for  x^  correspondingly  in 

X  pm  ,  .[t!°J  .In  fact,  it  is  easy  to  see  now  how  to  handle  the  more  general 

.  „  m+k,A  1-1,  l 
1=1 

situation  in  which  we  have  differential  operators  of  possibly  different  orders  on  the 
different  intervals  given  by  the  partition  AQ  ,  with  appropriate  side  conditions  at  the 
points  of  Aq  tying  the  pieces  together. 

Turning  now  to  the  eigenvalue  problem.  Corollary  1  of  Lemma  9.1  is  the  general  version 
promised  in  [2]  of  Lemma  3.1  there.  It  therefore  permits  the  following  generalization  of 
Theorem  3.1  there. 

Theorem  9.3  Let  T  =  NM  ^  be  the  compact  map  on  [  0  ,  1  ]  ,  1  <  p  *  00  ,  associated 

with  the  sufficiently  smooth  operators  M,  N,  and  0  of  (0.2).  Let  u  be  a  nonzero 

eigenvalue  of  T  with  corresponding  invariant  subspace  S  ,  and  let  J  be  a  matrix 

representation  for  T|^  .  Let  =  P^T  ,  where  is  the  projector  given  by  (9) 

associated  with  a  ssuper  projector  Q  of  order  (k,n).  Then,  for  all  small  |A|  ,  T^  has 

an  invariant  subspace  ,  and  T^  J  ^  has  a  matrix  representation  for  which 

A 

HJ  -  <  const  | A | 
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10.  bsuper  projectors  of  order  (k,k)  associated  with  point  evaluations.  We  now  loo 


in  more  detail  at  the  possible  ssuper  projectors  of  order  (k,k).  To  begin  with,  we  only 
consider  their  action  on  P2k  ,  and  this  we  can  describe  fully  by  specifying  their  action 
on  the  elements  of  some  basis  for  P2k  .  We  found  it  particularly  convenient  to  work  with 
the  basis  (p^)gk  1  consisting  of  the  Legendre  polynomials.  Then,  for  any  linear 
projector  Q  onto  Pk  , 


(1) 


QP .  =  P 

3  3 

k 

QP_  ■  =  1  a  ■  PL. 

k+3  r=1  gr  k-r 


g=0 , . . . ,k-1 


and  different  projectors  Q  correspond  to  different  matrices  (a^)  •  Further,  two  such 
projectors  agree  on  Pk+r  iff  the  corresponding  matrices  agree  in  rows  0,...,r-1  . 

Let 


f .  :  = 

3 


j-1 


j  <  k 


pj_1  "  ®Pj-1  '  >  k 


Then  (f .)  satisfies  (9.5a-b),  hence,  with  (9.5c),  Q  is  ssuper  of  order  (k,k)  iff 
,  k 

II  Pi-1rf1aj-k-1,rPk-r  =  0  for  1<k<^  * 

This  holds  iff  aj_k_i(r  =  0  for  i-1  =  k-r  and  i  <  k  <  j  <  2k+1-i  ,  i.e.,  iff 

a  =0  for  0  <  q  <  r-1 

V 

We  have  proved 


10.1  The  conditions 


(2) 


3-1 

Z  a  P  , 

r=1  3-1. r  k-r 

establish  a  one-to-one  correspondence  between  ssuper  projectors  Q  on  P2k 


Qf,  ,  .  =  o  with  f  :=  p  and  f,  :=  P 

k+3  -  k+1  k  -  k+j  k+j-1 


triangular  matrices  (a^J* 


j=2, . . . ,k , 
and  lower 


Of  course,  any  such  ssuper  projector  Q  on  P2k  can  be  extended  to  Infinitely  -nr.- 
such  on  C C— 1 ,11#  and  any  such  can  be  obtained  in  the  form  QP  ,  with  P  an  arbitral--.- 
linear  projector  on  C[-1,1]  onto  P2k  .  We  choose  to  ignore  this  aspect,  though.,  sin-i- 
the  property  of  beinq  ssuper  of  order  (k,k)  depends  only  on  the  action  on  P-,k  . 

-8  - 


Lemma  10.1  gives  rise  to  several  observations. 

The  first  interesting  basis  function,  viz.  »  is  simply  the  k-th  Legendre 

polynomial,  Pk  .  Hence,  if  we  think  now  of  Q  as  being  given  by  the  rule 

*  * 

(3)  Qf  E  and  q^f  =  qjf  ,  i=1,...,k 

*  * 

for  suitably  chosen  linear  functionals  q^ ,  ...,  q^  ,  then  we  must  have 

(4)  g*Pk  =  o  ,  i«i,...,k  . 

Now,  in  Example  1  (Collocation),  we  had 

* 

q.f  =  f ( >  ,  i=1,...,k 

and  so  (4)  is  satisfied  (for  n  =  k  in  Example  1)  since  then  <Pi>1  is  simply  the  sequence 
of  zeros  of  P.  .  In  Example  2  (Least  squares), 

q^f  -  /  P^_|f  ,  i=1,...,k  , 

and,  again,  (4)  is  satisfied  since  P^  is  orthogonal  to  =  span(Pi_1 .  Suppose  no* 
that,  in  an  attempt  to  bridge  the  gap  between  these  two  particular  ssuper  projectors,  we 
look  for  ssuper  projectors  for  which  some  of  the  interpolation  conditions  are  point 
evaluations,  say 

* 

qif  «  f(oi>  ,  i=1,...,r 

for  some  r  .  Then  we  conclude  from  (4)  that  {a^,  ...»  o^}  must  be  a  subset  of 

{p  ,  ...,p  }  :=  zeros  of  P,  .  This  leads  us  to  consider 
Ik  K 

Example  3.  Ssuper  projectors  using  point  evaluations  Let  (p^)^  be  tbe  zeros  of 
Pk  in  some  order.  Then,  for  r=0,...,k,  the  conditions 

(5)  O  f  6  P,  ,  f-QfiP  ,  Qf(p,)  =  f(p.)  for  i=r+1 , . . .  ,k 

r  k  r  r  r  i  i 

define  a  ssuper  projector  of  order  (k,k). 

This  provides  us  with  a  sequence  of  ssuper  projectors  of  order  (k,k),  with  Op  ,  i.»., 
interpolation  at  the  Gauss-Legendre  points,  at  one  end  and  ,  i.e.,  Least-squares 
approximation,  at  the  other,  and  so  demonstrates  a  perhaps  surprisinaly  simple  connection 


between  the  two 


We  no-;  verify  Example  1.  It  order  to  confirm  ^ )  dof  i  nes  1  linear  projector  2 

we  note  that  the  conditions  mentioned  are  equivalent  to  demand j nu  that 

*  * 

v»  f  e  P,  ,  q  .  Q  f  *  ti .  f  #  1  - 1  , .  • . ,  k 
r  k  i  r  i 

with 


(6) 


i 


m  pi-r 


e(p.  ) 

i 


i  :  1  ,  •  •  •  r 
i  -  r+ 1 ,  . . . ,  k 


Thus  it  suffices  to  show  that  the  matrix 
(7)  c 


is  invertible.  For  this,  assume  that  this  matrix  maps  a  -  [a.  y  '  to  0  ,  i.e., 

v 

* 

(8)  q^p  -  0  ,  i=1,...,k  with  p  :=  I 

Let  (w^)  be  the  woiqhc  vector  (known  to  be  strictly  positive)  for  the  corresnon.ii  n  i 
quadrature  rule 


-1 

Then  (6)  and  (8)  imply  that 

1 


l  w.f(p  .  )  ,  all  f.:P„ 

ri1  3  2k 


<  pp  .  .  *  E  W ,  p ( p . )  P .  ( P  .  )  ,  i  -  1  ,  • 

-i  ^  j.,  ^  3  i-i  i 


which  shows  that 


the  invertible  matrix  !p.  .  maps  the  vector  w  p(>.  »  *  to 

*  v  1  - 1  3  '  1  /  3 1  v  l  1 


0  and  consequently  p  not  only  vanishes  at  p 


r+ 1 

Thus  p  =  0  ,  and  so  a  -  0 
Note  that  the  invert ibi 1 i ty  of  (7)  just  proven  impli. 


p(p^)  =  0  for  j=1. 


(by  (ft)  and  (>*')  but  also 


i  nv»»rf  i  hi  l  j  tv  o 


(9) 


since  q. P  , 
i  3-1 


. P  .  , (P.  )  .  .  .  '  , 

3»1  1  '  1  ,  3-  r+  1 

P.  ,P.  ,  =  0  for  i  <*  r  <  i  . 

1-1  3-1 


To  verify  that  0r  is  ssupor  (a  Pact  not  int"i‘  li.U'»ly  *  1  us',  »vw  sh  '.%• 

that  2r  can  he  obtained  from  hy  a  suitable  nod  i  r  j  •*  i f  i-r;.  For  t'-is,  we  •>«».»  ]  t.o 

consi  ler  these  projectors  on  (on  v.'hi'h  th»»  r  pr  ''*•  *!  -rs  -if  >r  lor  !'•*,*  "i  *  u’"  1 

k  (  *  ♦  1 )  /•--  1  imens  innal  hynerpl  a  no )  *  T,.*t 


'1-?  )V 


\  ,  .  .  ,  ,  i* 


■  •  1  ,  .  .  .  ,  ♦  1 


h»*n*'  • 


-  1  > 


f$  =  pk+j-i  -  j/S-U  px-s  -  j=i . k+i 

[r]  sk-M.k 


for  some  matrix  (a  )  With  this  notation,  we  get,  for  r  =  0  , 

j  *  * s  3  * s—  I 


10] 

fk+j(Pi)  =  0  '  i=1'””k<  j=2,...,k+1 


and  so 


(10) 


X?1  .  ° 

p  (0  )  •  ••  p  (o  ) 

k-1  Mk  k-rp1; 

•  • 

•WV  —  pk+i(pi> 

•  • 

•  • 

•  • 

[0J  [0] 

•  • 

•  • 

= 

•  ■ 

•  • 

kl  kk  J 

WV  pk-k(pi) 

P  (p  )  •  •  •  p  ) 

k+k' pk ;  k+klp1IJ 

with  (a^  }  a  lower  triangular  matrix  (by  Lemma  10.1)  since  Q0  is  ssuper.  Now  write 
(101  in  terms  of  partitioned  matrices  as 


(10)' 


f  A.  0 

B  „  B 

‘c 

c 

11 

11  12 

_ 

11 

12 

A  ^  A_  _ 

B  B 

c 

c 

21  22J 

21  22 

*•  21 

22 

with  A11(  B11#  and  C1 1  all  of  order  k-r  .  Our  intent  is  to  replace  A^  by  0  and  to 
modify  A2?  correspondingly  in  such  a  way  that  the  equality  in  (10)  or  (10)’  is  preserved 
at  least  in  the  first  k-r  columns.  Explicitly, 


'21 


A2 1 B 1 1  +  A22B21 


(a2i  +  a22b2ibii  1 J  Bn 


(and  is  indeed  invertible  since  it  is  just  a  permutation  of  the  matrix  (9)).  Thus 


(11) 


with 


21 


An 

0 

B11 

rg 

00 

C1, 

C12 

A21 

0 

-B21 

B22 

C21 

^22 

:  = 

A21 

+ 

A22B21B11 

.  c22 

:  = 

A  2 1 B 1 2 

Now 


consider  the  linear  projector  Q  on  B2k+1  given  by 


(12) 


Qf . 


f  ,  i  <  k 
0  ,  i  >  k 


with 


(13) 


Pi-1 


,  i  <  k 


k+1  ■' 


and 


‘1 1 


J  3  ' 


k+i 


P,  .  -  Ha.  .P  .  ,  i  >  1 

k+i-1  _.  =  1  1-1,3  k-3 


'21 


Then  iva..)  is  lower  triangular,  hence  Q  is  a  ssuper  projector  of  order 


10.1  .  Further,  on  comparing  (10)— (10)*  with  (11),  we  see  that,  for  the  linear  f 


q .  of  (6 )  , 

l 


q.f  =  f,  (p.)  =  0  for  i=2,...,k+1,  j=r+1 , . 

^3  k+i  k+i  3 


while  the  fact  that  the  last  r  columns  of  (a_)  are  zero  implies  that 


i=2,...,k+1,  i.e.,  also 


"j  k+i 


In  addition,  trivially,  q  +  j  =  0  ,  all  j  .  We  conclude  that 


,  •»  k + 1  * 

ker  Ql  «  spani  f,  .1.  is  contained  in  hker  q.l  and  thus  must  eoua 

P„.  . .  k+3’l1  IIP, 


’ 2k+1  J  '  J  '  2k+1 

both  are  of  dimension  k+1  .  This  shows  that  Q  =  Qr  (on  P2k+1 


(14) 


(a[r]) 

13  ; 


11 


A21+fl22E21B1 1 


We  have  established,  in  particular,  that  2r  is  ssuper  of  order  (k,k1.  In 

comparing  again  (lO)-(IO)'  with  (11),  we  see  that  Qr  agrees  with  on  spa" 

1  P2k+1-r  ’  T^us  we  couldn't  tell  and  Qq  apart  on  p2k  * 

We  close  our  discussion  of  Example  3  with  the  following  four  observations. 

(i)  For  r  /  s,  differs  from  Q„  somewhere  on  P„.  ,  -  (w1 I  ' 

-  r  -  s  -  2  k  *2 -naxir,^- 

just  noted,  the  two  agree  with  Qp  ,  hence  with  each  other,  on  +  ^  ^ 

proof,  apply  both  sides  of  (10)  to  the  matrix 


j^Vl-i’  i  ,t;1 


■  k+  1-i  k 

Th»»n,  assuming  i  .at  the  Legendre  polynomials  have  all  h-»on  nirml. 


,  1  2 

1  *  r. 


i  i 


(» 


:  i n  1  that 


-11- 


(15) 

and,  in  particular, 

(16)  S 


[0] 

ii 


(a[°]  )  =  (  Z  w  P.  .  .  (p  )P.  .(p  ) 

1  ii  1  v  ,  s  k+i  s  k-i  s 

J  s=1 


^wsW^)Pk-i(PS1 

/-!pk+ipk-i  +  constkD2k(ViPk-i)  * 


This  shows  with  (14)  that  rank(aj^)  =  k-r  ,  all  r,  and  so  proves  our  assertion. 

(ii)  The  agreement  of  Qr  with  Qq  on  ^2k+1-r  not  restricte^  to  tlle  particular 

ssuper  projector  Qr  .  If  Q  is  any  ssuper  projector  of  order  (k,k)  which  enforces 
agreement  at  k-r  points,  then,  not  only  must  the  k-r  points  all  be  zeros  of  Pk  ,  say 
the  points  Pr+1»  •••»  ( in  some  suitable  ordering),  but  such  Q  then  necessarily 

agrees  with  Q0  on  P2k+1-r  *  For<  ^  I*emma  10*1,  5  satisfies  (2)  for  some  lower 
triangular  matrix.  The  matching  of  function  values  at  Pr+^  '  •• 
equality 


p^  then  forces  the 


a  0 

11 

•  • 

[Pk-1(pk>  —  Pk-l(pr+1>] 

P  (p  )  •  •  •  P  (p  ) 

k+1pk'  k+1lpr+1 

•  • 

•  • 

•  • 

•  • 

•  • 

•  • 

•  • 

A-r,1  ’**  Vr,k-t 

P  (£>)•••  P  ( p  ) 

k-r  Mk  k-r  r+1  J 

P„,  (P,  )•••  p,,  (P  , ) 

L  2k-r  k  2k-r  r+1 

In  the  terms  of  (10)',  this  reads 


I  1  k-r 


1  B11 


'1 1 


and  the  invertibility  of  used  earlier  now  proves  that  therefore 

i ■  i  k-r  _  _  .'  (0)1  k-r 

^  ij  Ji, j=1  ~  A11  "  ^  ij  i ,  j=1  ' 

(iii)  The  sequence  Q^,  ...»  Q^.i  connecting  QQ  to  Qk  constructed  in  Example  3 


depends  on  the  particular  order  p^. 


...»  p,  in  which  we  have  chosen  to  write  down  the 
k 


k  zeros  of  P^  .  _If  (Q^J  is  the  sequence  corresponding  to  the  ordering  p’,  .... 
then  Qr  =  (V  ^n  P2k+1  iff  the  two  sets  (Pf+1  ,  ...,  Pk>  and  {p^1#  •••»  P^; 


coincide.  Indeed,  from  (14),  on 

-1 


A21  +  A22B;i(B;i’ 


“1 


Now,  obviously#  A^  *  #  -  (15)  makes  this  quite  explicit,  hut  it  is  clear  anvway 

since  does  not  depend  on  the  order  in  which  we  write  down  the  interpolation  point - 


22 


2k+ 1 


-13- 


:=  B1  (?'  ) ~1  -  P  3-1 

2111  2111 


0  •  Let  .1  he  the  »*>•■* r.T  :t  -i  * 


Then,  on  partitioning  II  as  ft  is  in  (10)’,  wp  rir.  i 


?  ;;  +  p  .n’ 

nil  12  ji  '  *  21 


^21  “  B21B11B11 


=  B  .1  t  B  a  -  b  P  {  n  *  t  1 

2111  22  21  21111111  12  21 

’  (B22  -  B21BTlB12,::21 

How  note  that  ~  ^21B11  ^B12  *‘s  *'lle  lower  right  diagonal  block  obtained  by  block  ~ao<?s 

elimination  applied  to  B  ,  hence  is  invertible  (since  ft  in).  We  conduit  that  P  -  0 
iff  li^  -  0  ,  and  that  says  that  II  permutes  the  first  k-r  columns  of  ft  amor. a 
themselves . 

( iv)  Finally,  we  observed  earlier  that  the  collection  of  ill  ssuoer  projectors  of 
order  (k,k)  on  p2k+1  f°rms  d  linear  manifold  or  hyperplane  of  ii  mens  ion  k'k+U/2  .  \i** 
now  show  that  this  linear  manifold  is  spanned  by  the  particular  $ super  projectors  p 
introduced  here.  Precisely,  we  show  that  the  collection  og  ^n  ssnpor  projectors  of  order 
(k,k)  on  P2k+1  is  the  affine  hull  of  the  1  +  k(k+1)/2  particular  projectors 
Q0'  Qi  1  <  •  •  •  '21k.Q22'  •  •  •  '22H'Q33'  •  •  • 

V 

with  Qq  collocation  at  all  the  zeros  ( p ^ ^  of  ,  while,  for  1  <  r  <  a  <  k  , 

is  given  by  orthogonality  to  Pr  and  matching  of  function  values  at  the  k-r  points 

p  ,  P  „  ,  p  p,  .  Here,  (p.)  are  in  any  pa**t  i^-ilar  order;  i  i  fact,  in  the 

r  s-1  sm  .<  i 

definition  of  ,  thin  orl^r  could  even  change  with  r  (though  not  with  s  ).  To  nro;.- 

the  assertion,  it  is  sufficient  to  show  that  the  k(k+1)/2  1  inear  macs 
{o  -o  :  1  <r<;v»k  }  are  linearly  in  lonon-lont  (as  points  in  rh>*  linear  :pac.'  of  ill  linear 


naps  on  P?'+1  ^  and  for  this.  It  is  sufri«-i, 


i  nea v-  f met  iona  1  s 


n  t  t  '  exh  i  in  i  t 


e-  .  )  x  .  /  r' 

.  i  i  ] 


f  ■  r  i  >  r  »  r.  1 


since  this  insures  that  the  matrix  fy..(Q  -Q_)x.  1 

v  i]  rs  0  i]' 


is  upper  triangular  with  nonzero 


diagonal  entries  (using  the  ordering  1 1 , 1 2, . . . , Ik, 22, . . . , 2k, 33, . . . ,kk  ),  hence  invertible. 

(We  are  using  here  the  standard  argument  whereby  the  sequence  (y  )  in  a  linear  space  is 

linearly  independent  iff  there  exists  a  corresponding  sequence  (vr)  ° *  linear 

functionals  on  that  space  for  which  the  matrix  (v^y^)  *s  invertible.)  First,  pick  x^j  = 

p2k+i-i<  •  Then,  since  Qrg  forces  agreement  at  k-r  points,  it  agrees,  by  (ii), 

with  Q0  on  P2k+1-r  and  so  w i j (Qrs”®0 *xij  =  0  for  2k+1_i  <  2k+1-r  ,  i.e.,  for  i  >  r 

no  matter  how  we  pick  y^  .  Further,  pick  y_:  ^  <  all  i/j  •  Then,  as  both 

Qn  and  Q__  match  the  value  at  p.  when  r  <  j  ^  s  ,  we  conclude  that  y.  .(Q  -Q n)xf.  = 
u  rs  j  j  j  rs  u  i  j 

0  also  for  i  =  r  and  j  >  s  (>r)  .  Finally,  we  claim  that  y^ j”^0 *xi j  ^  0  '  For' 

otherwise,  2ijP2k+1-i  would  agree  with  p2k+1-i  at  Pi'  •••»  Pk  as  wel1  as  at  the 

linear  functionals  f  1 — *  / P^.  ,  r=1,...i-1,  i.e.,  would  agree  with  Qi_1  i_1  at 

P2k+1-i  an<^  k'1^s  woul‘i  contradict  (i). 

Finally,  up  to  this  point,  this  section  has  been  concerned  with  ssuper  projectors  of 
order  (k,k).  But  we  think  it  worth  recording  a  version  of  Lemma  10.1  for  ssuper 
projectors  Q  of  order  (k,n)  ,  1  <  n  <  k  ,  along  with  a  corresponding  corollary  concerning 

t  * ik 

the  k  linear  functionals  associated  with  Q  . 

Lemma  10.2  The  conditions 

(17)  2fk+j  =  0  with  fk+j  :=  Vj.,  -  C"3"1  aj-1,rPk-r  *  ^ . N 

establish  a  one-to-one  correspondence  betwen  ssuper  projectors  Q  of  order  (k,n)  on  pk+N 
,  N  >  n,  and  lower  trapezoidal  matrices 

(18)  (a  )N  }  \  ,  with  a  =0  for  k-n+q  <  r  <  k  . 

v  qr'q=0  r=  1  -  qr  - 

.  *  ik 

In  terms  of  linear  functionals  ( J  ^  associated  with  Q  via  (3),  the  ssuper 
projector  criterion  (9.4)  may  be  expressed  as 


p.  i 


rv- 

r  -  1 


an  i  lenmu  10.2  rray  be  restate*!  d;-; 


Corollary  The  linear  functionals  n  ^  define 
^  PkVi 


k*k 


'Nxk 


with  b . 


ij 


q  .  P ,  ,  .  ,  0  <  i  <  k  /  d  •  i 

3  k-1-l  - 


(whose  transpose  decribes  the  action  of  the  functionals 

•••'  P0'  Pk . Pk+N-1>  of  Pk  +  M  >  sdti-?-fX 

CS*1  =  .  A  , 

■jherfl  the  matrix  A  =  fa  .  has  the  lower  traneroi  >V 

- - -  —  v  IV  1=0  i=1  - - - - - - 


functionals  define  distinct  ssuper  projectors  precisely 
corresponding  matrices  A  are  distinct. 


Proof.  Q  is  a  projector  onto 
(k,n)  on  iff  the  matrix  A 

element  of  the  assertion  C  -  AB  = 


Pk  iff  B  is  nonsingul ir. 
of  (18)  is  connected  with 

* 

0  is  found  by  apply inu  ..i 


-1  - 


! 
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